Spin-polarized Josephson and quasiparticle currents in superconducting spin-filter tunnel junctions 
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We present a theoretical study of the effect of spin-filtering on the Josephson and dissipative quasiparticle cur- 
rents in a superconducting tunnel junction. By combining the quasiclassical Green's functions and the tunneling 
Hamiltonian method we describe the transport properties of a generic junction consisting of two superconduct- 
ing leads with an effective exchange field h separated by a spin-filter insulating barrier. We show that besides the 
tunneling of Cooper pairs with total spin-projection S z = there is another contribution to the Josephson current 
due to triplet Cooper pairs with total spin-projection S z ^ 0. The latter is finite and not affected by the spin-filter 
effect provided that the fields h and the magnetization of the barrier are non-collinear . We also determine the 
quasiparticle current for a symmetric junction and show that the differential conductance may exhibit peaks at 
different values of the voltage depending on the polarization of the spin-filter, and the relative angle between the 
exchange fields and the magnetization of the barrier. Our findings provide a plausible explanation for existing 
experiments on Josephson junctions with magnetic barriers, predict new effects and show how spin-polarized 
supercurrents in hybrid structures can be created. 

PACS numbers: 74.50.+r, 72.25.-b, 74.78.Fk 



Introduction The prediction of long-range triplet supercon- 
ducting correlations in superconductor-ferromagnet (S-F) hy- 
brid structures!^ has led to intense experimental activity in 
the last years 3 - - — . These experiments have shown that a fi- 
nite Josephson current can flow between two superconductors 
connected by a ferromagnetic layer whose thickness far ex- 
ceeds the expected penetration length of singlet pairs. The 
Josephson current measured in these experiments is attributed 
to the flow of Cooper pairs in a triplet state. According to 
the theory, the appearance of triplet correlations occurs only 
in the presence of a magnetic inhomogeneity located in the 
vicinity of the SF interfac o'i 8 ^ 11 ' 34 . The inhomogeneity can 
be either artificially created^ or can be an intrinsic property 
of the material, as for example the domain structure of usual 
ferromagnets 6 - or the spiral-like magnetization in certain rare- 
earth metals^i. 

The Josephson triplet current is nothing but a dissipation- 
less spin-polarized current and therefore its control would 
be of great advantage in the field of spintronics 12 . Impor- 
tant building blocks of spintronic circuits are magnetic in- 
sulating barriers with spin-dependent transmission, so called 
spin-filters (Is/), which have been studied in several ex- 
periments using for example europium chalcogenide tunnel 
barriers^ " 15 ! 30 The question naturally arises whether one can 
use these spin-filter tunneling junctions to control and eventu- 
ally to create a triplet Josephson current. We will address this 
question in the present letter. 

In spite of several studies of the transport properties of 
spin-filter tunneling barriers, the Josephson effect has only 
recently been explored through a S-I s /-S^. The tunnel bar- 
rier used was a GdN film that reduced the value of criti- 
cal current I c compared to a non-magnetic barrier. Beside 



the large reduction of I c the authors of RefJ£ also observed 
that the I C (T) curve deviates at low temperature from the ex- 
pected tunneling behavior 20 . Teoretically, the effect of spin- 
dependent transmission on the Josephson current was first 
considered by Kulik 2 -! and Bulaevskii et al^ 2 . on the basis of 
the tunneling Hamiltonian. It was demonstrated that spin- 
selective tunneling always leads to a reduction of the crit- 
ical current with respect to its value in the spin indepen- 
dent case or even to the change of sign of the critical cur- 
rent. Later on it was shown that the magnetic barrier in 
a Ij/-S structure induces an effective exchange field in the 
superconducto r 16 ! 17 . Other theoretical works have addressed 
the Josephson effect through spin-active barriers in ballis- 
tic systems^ 3 -"— and through ideally ballistic superconductor- 
ferromagnetic insulator-superconductor junction s 26 ' 27 . Also 
the spin-polarized current through S-N-F junction has been 
studied in Refs i 28 ' 29 . However, none of these works presented 
a comprehensive theoretical study of the Josephson effect by 
taking into account both the spin-filter effect and the presence 
of the exchange field in the superconducting electrodes. Nor 
has the interplay between spin-filtering and triplet supercur- 
rents been investigated. 

The aim of the present letter is to provide a complete de- 
scription of the transport properties of Josephson junctions 
with spin-filters. For that sake we introduce a simple model 
which allows us on the one hand to derive simple and useful 
expressions for the dc Josephson and quasiparticle currents 
in a S-I s f-S and on the other hand to predict the conditions 
under which the creation of a spin-polarized super-current is 
possible. Our model considers the spin-filter effect of the 
Is/ barrier and a finite exchange field in the electrodes. We 
show that the contribution to the current from Cooper pairs 
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in the singlet and triplet states with zero spin projection van- 
ishes in the case of a fully spin-polarized barrier. However, 
the contribution to the Josephson current from tunneling of 
Cooper pairs in a triplet state with non-vanishing spin pro- 
jection is independent of the strength of the spin-filter. The 
latter contribution is finite provided that the exchange fields 
in the electrodes and the spin quantization axis of the barrier 
are non-collinear. This remarkable result explains how spin- 
polarized currents can be created and controlled by means of 
spin-filter barriers. We also calculate the differential tunnel- 
ing conductance of the S-I s /-S junction, and analyze how the 
Zeeman-splitting peaks depend on both the spin-filter parame- 
ter and the exchange field in the electrodes. Our model allows 
a quantitative description of existing transport experiments on 
S-T/-S junctions^"— and gives a possible explanation for the 
temperature dependence of the critical Josephson current ob- 
served in Ref.— , We finally discuss the applicability of our 
model to real systems. 

The Model We consider a tunnel junction between two su- 
perconductors (see inset of Fig. 1). The tunneling barrier, 
black area in Fig. 1, is a spin-filter. The grey regions close to 
the barrier are thin ferromagnetic layers with a finite exchange 
field acting on the spin of the conducting electrons. The di- 
rection of these fields is arbitrary. We assume for simplicity 
that the thickness of the superconductors is smaller than the 
coherence length. In this case one can average the equations 
for the Green functions over the thickness and get a uniform 
superconductor with built-in exchange field 3 ^. Under these 
assumptions the system is described by a generic Hamiltonian 
which is homogeneous in space: 

H = H R + H L + H T (1) 

where H R ^ describes the left and right electrodes consisting 
of a BCS superconductor with an intrinsic exchange field. For 
example for the left electrode it reads 

h l = E a ls IZA' ~ ( /! L n L-6).«'] flfa+E ( A L a 'ki a -ki + hx ) 

k,s,s' k 

(2) 

where a{a i ) are the annihilation (creation) operator of a parti- 
cle with momentum k and spin s, ^ is the quasiparticle energy, 
A is the superconducting gap, a = (<7i , 02, 03) is a vector with 
the Pauli matrices, Iil the amplitude of the effective exchange 
field and n a unit vector pointing in its direction. The Hj term 
in Eq. (Q3 describes the spin-selective tunneling through the 
spin-filter and is given by 

H T = £ (7& + U6 Z ) SS , a\b s , + h.c. (3) 

s,s' 

where a and b are the field operators in the left and right elec- 
trodes respectively. 7 and l( are the spin independent and 
spin dependent tunneling matrix elements. We neglect their 
momentum dependence. The tunneling amplitude for spin up 
(down) is then given by T^m —7 ±11. We assume that the 
origin of the different tunneling amplitudes is the conduction- 
band splitting in the ferromagnetic insulating barrier which 
leads to different tunnel barrier heights for spin-up and spin- 
down electron s 15 ' 30 . 



In order to calculate the current through the junction it 
is convenient to introduce the quasiclassical Green functions 
Sr(l) f° r the left and right electrodes. An expression for the 
current in terms of g R {j\ can be obtained straightforwardly 
from the equations of motions for the Green functions after 
integration over the quasiparticle energy. In the lowest order 
in tunneling the current is given by 

I= 32^ + 1^) 1^^^} (4) 

where R N = l/[47t{eN{0)) 2 {T^ + Tj 2 )] is the resistance of the 
barrier in the normal state, N(0) is the normal density of states 
at the Fermi level, the symbol T denotes 8x8 matrices in the 
Gor'kov-Nambu (t,) - Spin (a,) - Keldysh space, f = 7% (g> 
Ob + ® 63, a, j3 are the angles between the exchange field 
of the L and R electrode with respect to the z-axis (see inset 
in Fig. 1), and g a is the bulk Green's function which can be 
obtained by solving the quasiclassical equations. The matrix 
g La (and in analogy g R p) can be written as g La = Ra-guaMi, 
where gia is the known solution for the case of an exchange 
field along the z— axis, andR a = cos(a/2) + 61 sin(a/2) 
Results We first proceed to determine the Josephson critical 
current through the spin-filter. We assume that A^ = A R = A 
and hi = Iir = h and that the exchange field in the left (right) 
electrode forms an angle a (j3) with the magnetization of the 
Is/ barrier which points in z-direction. From Eq. (0]i we find 
the Josephson current Ij = I c sin 9, where <p is the phase dif- 
ference between the superconductors and the critical current 
I c is given by the general expression: 

eR N I c = 2itT £ {r [f 2 + / 2 cosacosj3] + / 2 sin a sin j3 } , 

fl)„>0 

(5) 

where r = IT^^T 2 + T 2 ) is a parameter describing the 
efficiency of the spin-filtering (r = denotes full polariza- 
tion while r = 1 a non-magnetic barrier™) . f s t t \ = (/+ ± 
f-)/2 are the anomalous Green's functions where f± — 
A/ \J {o) n ± ih) 2 ±A 2 and (O n the Matsubara frequency. The 
amplitude of the singlet component is determined by f s 
whereas the amplitude of the triplet component is given by 
/ f . Eq. (O is one of the main results of our work. If h = 
it reproduces the expression presented in Refs i 21 ' 22 which is 
the well-known Ambegaokar-Baratoff (AB) formula for the 
critical current^ , multiplied by a factor r < 1 . 

In the case of a fully spin-polarized barrier (r = 0), i. e. if 
either 7| or 7+ is zero, Eq.© shows that the singlet Cooper 
pairs do not contribute to the Josephson current. The contri- 
bution to the current is only due to the second term on the 
r.h.s. which is independent of r and proportional to the am- 
plitude of the triplet component f t . This term does not vanish 
provided that neither a nor /3 are equal to or %. This im- 
portant result shows that even though in the electrodes only 
the triplet component with (locally) zero spin projection ex- 
ists, the non-collinearity between h and the magnetization of 
the barrier induces a coupling between them and leads to the 
creation of a spin polarized supercurrent. In our model the 
parameters r and h are independent. However, for a ferromag- 
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FIG. 1. (Color online) The temperature dependence of the critical 
current for different values of h and r. We assume that a = j8 = 0. 
Inset: The structure described by our model Hamiltonian Eq. (0. 
The black region represents the spin-filter barrier while the grey re- 
gions are layers with a finite exchange field pointing in arbitrary 
direction, dp is the thickness of these layers and we have defined 
I c = 2I c eRj^/(Aon), where Ao is the value of the order parameter at 
T = and h = 0. 



netic insulator/ superconductor system they might be related 
to each other—. 

We assume next that the exchange fields in the left and 
right electrodes are parallel to the magnetization of the bar- 
rier (a = j3 = 0) and compute the temperature dependence of 
the critical current using Eq. (0. In panel FigQ] we show this 
dependence for different sets of parameters (h,r). Through- 
out this article the order parameter A(T) is determined self- 
consistently and the temperature in the figures is normalized 
with respect to the critical temperature which depends on h. 
The I C (T) curve was measured in Ref- 8 - for a Josephson junc- 
tion with a spin-filter as tunneling barrier. If we assume, as the 
authors of Ref.— did, a finite spin-filtering effect (r < 1) but 
neglect the exchange field in the superconductor (h = 0) we 
obtain the dashed curve in FigQ] which is nothing but the AB 
curve multiplied by a pre-factor r ps 0.27. If we assume now a 
finite value of the effective exchange field in the S layers (dot- 
dashed curve in FigJTJ for the same value of r one obtains that 
critical current is for all values of temperatures smaller than 
the AB curve. If we now keep the same value for the finite 
exchange field but slightly change the value of r ( solid line in 
Fig|T), one can see that for lower temperatures the I C {T) curve 
exceeds that of the AB curve. This behavior, which is in qual- 
itative agreement with the results of Ref^, shows that the 
interplay between h and r is crucial to understand the trans- 
port properties of the junction. We cannot conclusively say, 
though, that the experiment can be fully explained by these 
results. Indeed, measurements of the tunneling conductance 
injunctions with GdN barriers suggest a finite exchange field 
inside the S electrodes 19 . However, GdN barriers may also ex- 
hibit a complicated temperature dependent magnetic domain 
structure that could also modify the I C {T) behavior—. This 
hypothetical effect is beyond the scope of the present work. 

Let us now assume that the exchange fields in the S lay- 
ers and the magnetization of the l s t barrier are non-collinear 
(we set a = j3 = n/2 in order to maximize the contribution of 
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FIG. 2. (Color online) The temperature dependence of I c for different 
values of r and a = 7t/2 (a); and for different values of a and r = 0.1 
(b). In both panels h = 0.567Ao 

the triplet supercurrent). In Fig. [2^ we show the temperature 
dependence of the critical current for different values of the 
spin-filter parameter r corresponding to highly polarized bar- 
riers. For large values of r the critical current is positive for 
all temperatures (0-junction). However if r is small enough 
the second term in the r.h.s. of © start to dominate and at 
certain interval of temperature I t < (n -junction), i.e. our 
model predicts a zero-7T transition for large enough spin-filter 
efficiency. Thus, it is more likely to observe the 0-n transition 
in systems containing europium chalcogenide tunnel barriers 
with a almost 100% spin polarization 30 than using GdN films 
with a spin-filter efficiency of around 75%' - . Note that in the 
fully polarized limit r = the critical current is negative for 
all temperatures. In Fig. [2J5 we show the I C {T) dependence 
for r = 0.1 and different values of a. Negative values of the 
current appear if a is close to n/2. The origin of the n junc- 
tion behavior described here is different from the one studied 
in Ref.—. The n -junction behavior shown in Fig. 2 is caused 
by the non-collinearity of the exchange fields and the mag- 
netization of the I s f, i.e. it is determined by the second term 
in the r.h.s of Eq. ©. In contrast, in Ref. 22 there is no such 
term and the n junction behavior was obtained by assuming 
that 1 < U (i.e. by choosing r < 0). For completeness we 
note that the Josephson current in metallic multilayered SFFS 
junctions also depends on the angle of mutual magnetization 
orientations in different F layers. This problem (without spin- 
filter barriers) was studied in numerous papers on the basis of 
Usadel, Eilenberger or Bogoliubov-de Gennes equations (see 
for exampl e) 8 ' 9 ' 35- ^ and references in the review articles^ and 

m. 

Let us now calculate the quasiparticle current I qp from 
Eq.©. For the normalized current j qp = I [/p (V)/In(V) 
Un(V) — V /Rn is the current through the junction in the nor- 
mal state) we get 

jap = J deF v Y aP (e,h,V) , (6) 

where Y a p(e,h,V) is the spectral conductance and Fy = 
0.5{tanh[(e + eV/2)/2T] - tanh[(e - eV/2)/2T}}. We 
present here the expression for a symmetric junction, i.e. i.e. 
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FIG. 3. (Color online) The zero- temperature normalized differential 
conductance G qp = Rj^dI qp /dV for h = 0.4Aq and r = 0,0.5, 1. Left 
panel: a = j3 = Jr/4, right panel: a = 0, jS = n. By calculating the 
curves we have added a small r\ = O.OlAo damping factor. Note that 
the curves are shifted vertically for clarity 



v, = V/ and a = j3, although similar expressions hold for ar- 
bitrary angles a and j3 . It reads, 

Y aa = Vq+Vq- + v 3+ v 3 _ - (1 - r)v 3+ v 3 _ sin 2 a . (7) 

We have denned the DOS v , 3 (e) = [v(fi + h) ± v(e - h)]/2 
with v(e) = e/Ve 2 -A 2 and v ± = v (e±eV/2) . 

The left panel of Fig. 4 shows the voltage dependence of the 
normalized differential conductance G qp for zero-temperature 
in the symmetric case a = j3 = n/A. In the absence of spin- 
filter effect (r = 1) the differential conductance is an even 
function of V, showing a peak at eV = 2A and no signature 
of the exchange splitting^ 2 - (dash-dotted line in the left panel 
of Fig. 4). However for r < 1 and < a < n two addi- 
tional peaks appear at eV = 2(A±/i). Notice that the height 
of these peaks increases with decreasing r. Thus, by measur- 
ing the differential conductance one can extract information 
about the model parameters a, r and h. From Eqs. (|6]|7]i one 
can also show that in order to obtain an asymmetric G qp (V) 
dependence one should set ^ Iir as discussed RefJ^. In 
the antiparallel case, where a = 0, j3 = % (right panel of Fig. 



4) and for r < 1 the differential conductance has peaks at 
eV = 2(A± h) (but not at eV = 2A). These peaks have the 
same size for r = 1 . However, by decreasing r towards zero, 
the difference between the peak sizes increases. In the fully 
spin-polarized case (r = 0) one of these peaks vanishes. 

Discussions and Conclusions Our model can describe dif- 
ferent systems. First, the model applies to junctions made of 
two magnetic superconductors separated by a spin-filter bar- 
rier 7 s f. But it can also describe a S-F-I s f-F-S junction with 
the width of the F-S electrodes smaller than the characteris- 
tic length over which the Green functions vary. We have also 
verified that our results are qualitatively valid for long S elec- 
trodes. These results will be discussed in more detail else- 
where. Finally, our model can also describe a simple S-I s f-S 
assuming that the l s f barrier induces an effective exchange 
field in the superconductor over distances of the order of the 
superconducting coherence length as predicted in Ref.— . 

In conclusion, by combining the quasiclassical Green func- 
tions and the tunneling Hamiltonian approach we have studied 
the effect of spin filtering on the Josephson and quasiparticle 
current in tunneling junctions. We have shown that for fully 
polarized barriers the singlet component does not contribute to 
the supercurrent/j. However, if the direction of the exchange 
field h in both electrodes is not parallel to the quantization 
axis of the barrier, a non-zero Ij current is observed due to 
the triplet component. In this case the current is 100% spin 
polarized. We have also calculated the differential conduc- 
tance and shown its dependence on the spin-filter parameter 
r and the misalignment angle. By measuring the differential 
conductance one can extract information about the magnetic 
structure of the spin-filter junction. Our findings are relevant 
for the creation, control and manipulation of spin polarized su- 
percurrents as well as for the characterization of S-F/-S junc- 
tions. 
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